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ABSTRACT 

At the final stage of terrestrial planet formation, known as the giant impact 
stage, a few tens of Mars-sized protoplanets collide with one another to form 
terrestrial planets. Almost all previous studies on the orbital and accretional 
evolution of protoplanets in this stage have been based on the assumption of per- 
fect accretion, where two colliding protoplanets always merge. However, recent 
impact simulations have shown that collisions among protoplanets are not always 
merging events, that is, two colliding protoplanets sometimes move apart after 
the collision (hit-and-run collision). As a first step towards studying the effects 
of such imperfect accretion of protoplanets on terrestrial planet formation, we 
investigated the merging criteria for collisions of rocky protoplanets. Using the 
smoothed particle hydrodynamic (SPH) method, we performed more than 1000 
simulations of giant impacts with various parameter sets, such as the mass ratio 
of protoplanets, 7, the total mass of two protoplanets, Mt, the impact angle, 6, 
and the impact velocity, fimp- We investigated the critical impact velocity, fcr, 
at the transition between merging and hit-and-run collisions. We found that the 
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normalized critical impact velocity, Vcv/vesc, depends on 7 and 6, but does not 
depend on Mt, where fesc is the two-body escape velocity. We derived a simple 
formula for fcr/^^esc as a function of 7 and 9 (Eq. flTB]) ). and applied it to the 
giant impact events obtained by A'^-body calculations in the previous studies. We 
found that 40% of these events should not be merging events. 

Subject headings: accretion, accretion disk — planets and satellites: formation 
— solar system: formation 



Introduction 



Planets are formed in a disk around a star called a protoplanetary disk, which is com- 
posed of gas and dust. Terrestrial planets are formed mainly from the dust component. 
Their formation process can be divided into three stages. The first stage is the formation 
of a large num ber of kilometer sized bodies called planetesim als by accretion among dust 
particles (e.g., iGoldreich fc WardI Il973l : lYoudin fc Shul |2002| ). In the second stage, these 
planetesimals collide to produce a fe w tens of Mars-sized objects called protoplanets (e.g., 
Wetherilll Il985l : iKokubo fc Idalll998r). The final stage is the formation of terrestrial plan- 
ets from protoplanets (e.g.. Ichambers fc Wetherill 1998 : Agnor et al. 1999 ). The collisions 
among protoplanets are referred to as giant impacts, and thus this final stage is known as 
the giant impact stage. 

Giant impacts have a large influence on the various features such as the number of 
terre strial planets formed, their mass and spin state (e.g., lAgnor et al.lll999l : iKokubo et al. 
20061 ). Giant impacts are highly en ergetic events , and are responsible for the creation of 



large satellites, like th e Moon fe.g.. ICanupl l2004al ) and planets with extremely large cores 
such as Mercury (e.g., Benz et al. 20071 ). Mo reover, giant impacts a re closely related to the 



thermal state such as a magm a ocean (e.g.. 



terrestrial planet atmospheres (iGenda fc Abe 



20051 ). 



bnks fc MeloshI Il992l ) , and the origins of the 



A large number of simulations of giant impacts have been devoted to the specific gi 
ant impact events related to the origin of the Moon or Mercury. However, since it is re- 
cently believed that multiple giant impacts are comnion dur i ng the last stage of terrestrial 



planet form ation, several studies (lAgnor &: Asphaua l2004al : lAsphaugl l2009l : iMarcus et al. 



2009 



ters. 



2010h have invest i gated t he giant impact simulations under various impact parame- 



Agnor fc Asphaugl (l2004al ) were the first to show that collisions of protoplanets during 



the giant impact stage are not always merging events, that is, two colliding protoplanets 
sometimes move apart after the collision. They called such a collision a hit-and-run collision. 
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Except for lKokubo fc Gendal (120101 ). all the previous studies on the orbital and accretional 
evolution of protoplanets during the giant impact stage have been based on the assumption 
of perfect accretion, where two colliding protoplanets alw ays merge. However, the hit-and- 
run collisions demonstrated by lAgnor &: Asphaugl (l2004al ) may have an important influence 
on many of the physical characteristics of terrestrial planets. 

In order to investigate the effects of such imperfect accretion of protoplanets on terres- 
trial planet format i on, the merging criteria during protoplanet collisions must be clarified. 
Agnor fc Asphaugl (j2004al ) performed 48 simulations of the collisions between same-sized 
protoplanets with masses of O.IM®, where is the Earth mass. They found that hit-and- 

for ^ = 30°, and Vi^p > 1.2vesc for 9 = 45° or 



run collisions occurred when Vimp > 1.5f, 



60°, where Vi^p is the impact velocity, fesc is the two-body escape velocity (equation [5]), 
and 9 is the impact angle. Since they varied fimp in steps of O.lfgsc for low-velocity col- 
lisions, the transition between the merging and hit-and-run collisions was estimated to be 



^imp 



1.4 - 1.5v,,r for 9 = 30°, and n 



imp 



1.1 - 1.2^;,,, for 9 = 45° or 60° 



Subsequently, I Agnor fc Asphaugl (l2004bl ) investigated collisions between diffe rent-sized 
protop lanets with mass ratios of 1:2 and 1:10. Their results were presented in lAsphaug 
(120091 ). The transition was able to be estimated as follows. In the case of a 1:2 mass ratio. 



the transition occurs at v- 



imp 



1.4 - 1.5t;psc for 9 = 30° and v 



60°, whereas in a case of 1:10 mass ratio, v 
for 9 = 45°, and Vi 



imp 



1.5-2.0t;e 



imp 

for 9 



l.l-1.2t;, 
= 30°, n 



imp 



for 9 = 45° or 
= 1.2- 1.3t;f,sc 



imp 



1.0 - l.lv,,r for 9 = 60°. 



Marcus et al.l ( l2009l ) and iMarcus et al.l ( l2010l ) investigated collisions of the rocky and 



icy super-Earths (up to lOM®), respectively. They focused on the stripping of the rocky or 
icy mantle resulting from a high- velocity impact. Although they did not perform the various 
simulations with low-impact velocities, we were able to find that the transition occurs at 



roughly v 



imp 



1.0 — 1.5fesc from their figures. 



As described above, the merging criteria has already been roughly determined for certain 
discrete values of the impact parameters. However, in order to carry out N-body orbital 
calculations with the merging criteria for the giant impact stage, a simple formula describing 
the dependence of the merging criteria on the impact parameters is required. To achieve this, 
it is necessary to determine the merging criteria over a wide range of impact parameters. 
In the present study, we performed more than 1000 simulations of giant impacts for various 
impact parameter sets using the smoothed particle hydrodynamic (SPH) method in order 
to formulate the merging criteria. 

In Section |2l we present the SPH code and initial conditions used in our giant impact 
simulations. In Section |3l we show the colhsion outcomes, and investigate the transition 
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between merging and hit-and-run collisions. We also perform a resolution test on the simu- 
lations. In Section |U we derive the merging criteria as a function of the impact parameters, 
and compare the results with those of previous studies. Using the derived criteria, we then 
discuss the merging probability of protoplanets during the giant impact stage in Section |3 



Calculation Method 



2.1. Numerical Code 



In order to p erform impact simulations for protoplanets, we used the SPH method (e.g., 
Monaghanlll992l ) . which is a flexible Lagrangian method of solving hydrodynamic equations, 
and has been widely used in previous giant impact simulations. The SPH method ca n easily 
deal w ith large deformations and shock waves. Our numerical code is based on ICanup 
(l2004bl ): here, we briefly describe its essential points. 



The equation of the motion for the i-th SPH particle is given by 



~dt 



all 



neighbor 



where Vj is the velocity of the i-th SPH particle, t is the time, and J^ij and Qij are the 
pressure gradient and mutual gravity terms between the z-th and j-th particles, respectively. 
Several forms have been used for the pressure gradient term, none of which appear to be 
clearly superior to the others. In this paper, we use the following symmetric expression. 



4 + 4 + H 



ViW{rij,h 



(2) 



where rrij, Pj, and pj are the mass, pressure, and density of the j-th particle, respectively, 
Hjj is the artificial viscosity, W is the kernel function, r^j is the distance between the z-th 
and j-th particles, and hij is the average smoothing length of the i-th and j-th particles. For 
the artificial viscosity Hjj, we use a Von Neur nann-Richt r nyer-t ype viscosity with parameters 
of Q^vis = 1-5 and /3vis = 3.0, as described in iMonaghanl (Il992l ). For the kernel function W, 
we use the spherically symmetric spline kernel function proposed by iMonaghan fc Lattanzio 
Jl985h : 



1 



W{r, h) 



3/r ■ 



2 _|_ 3/rN3 



K2- 



0<f <1, 

1 < 71 < 2, 

2 < f . 



(3) 



This function satisfies J W^(|r|, h)dY = 1, and has a zero value when r > 2h. In our code, h 
is variable for each particle and with time, and determined to satisfy the condition that the 
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number of neighboring particles (A^nci) within 2h is almost constant, A^nci = 64±2. We used 
a maximum value of the smoothing length (/imax) to save computational cost when searching 
for neighboring particles and calculating the pressure gradient term. The value of /imax is 
determined from 

^(2/imax)', (4) 



Pmin 3 

where pmin is the minimum density, which is set to 5 x 10^'^ kg/m^ in our simulations. 
The mutual gravity term in equation ([1]) can be written as 



gY.^"-^. (5) 



where G is the gravitational constant, and rhj is the effective mass of the j-th particle toward 
the i-th. particle defined by 



rhj = I Arcr mjW{r, hj)dr. (6) 





This equation gives rfij = rrij when rij > 2hj. 

The mutual gravity term between all SPH particles was directly comp uted using a 



speci al-purpose computer for gravitational N-body systems named GRAPE-6A (IFukushige et al 



20051 ). The GRAPE-6A can search for and produce lists of neighboring particles while si- 
multaneously calculating their mutual gravity. This list of neighbors is used to compute 
the pressure gradient term in the equation of motion and the time derivative of the internal 
energy. In a simulation over a period of 10^ sec with 20,000 SPH particles, the typical CPU 
time was about 4 hours. Thus, the GRAPE-6A allowed us to systematically explore a wide 
range of impact parameters. Time integration was performed using a PEC (predict, evaluate. 



and correct) scheme with variable time steps (e.g.. lSerna et al.lll996l ). which is second-order 
accurate in time. 



2.2. Pre-impact Protoplanets 



Here we describe the modeling method for the pre-impact protoplanets. All the proto- 
planets are assumed to be differentiated, with a 30% iron core and 70% silica.te ma ntle by 
mass. In our SPH simulations, we used the Tillotson equation of state (lTillotsonlll962n. which 



has been widely applied to giant impact simulations invo l ving shock waves (e.g . 



Benz et al. 



19871 : ICanup fc Asphaud I2OOII : lAgnor fc Asphaud l2004al : lAsphaue et al.ll2006f ). The Tillot- 



son equation of state contains ten material parameters, and the pressure is expressed as a 
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function of the density and the specific internal energy, which is convenient for treating fluid 
dynamics. We used the parameter se ts of granite fo r the silicate mantle and iron for the iron 



core, which are listed on page 234 of iMeloshl ( 119891 ) . 



All SPH particles in the protoplanets was set to have the same mass, and the total 
number of particles used for impact simulations was fixed at 20,000. For example, in the 
case of a collision of protoplanets with a mass ratio of 1:9, the smaller protoplanet consisted 
of 2,000 particles, and the larger one 18,000 particles. To model the pre-impact protoplanets, 
we placed the SPH particles in a 3D lattice (face-centered cubic) with iron particles on the 
inside, and rocky particles outside. The internal energy of the SPH particles was set to 
1.0 X 10^ J/kg. 

Beginning with this configuration, we calculated vibrations of the protoplanet until the 
particle velocities become slower than 100 m/s, which is much less than the impact velocity 
(order of km/s). After this operation, we used these relaxed objects as the protoplanets for 
impact simulation. As a first step, the protoplanets were assumed to have no spin. 

In order to set impact parameters such as the impact velocity (fimp) and impact an- 
gle (6*), we need to determine the radius of the pre-impact protoplanets. Since the surface 
boundary of the protoplanets described by the SPH particles is obscure owing to the smooth- 
ing length (h), we determined the radius of the pre-impact protoplanet (Rp) by the following 
equation. 



2.3. Initial Conditions for Collisions 

We prepared more than 1000 sets of initial conditions for the giant impact simulations. 
The parameters used were the mass ratio of the protoplanets (7 = Mi/Mt, where Mt and Mj 
are the mass of the target and impactor, respectively), the total mass of the two protoplanets 
(Mt = Mi + Mt), the impact angle (6*), and the impact velocity (fimp)- We systematically 
varied the mass ratio as 7 = 1, 2/3, 1/2, 1/3, 1/4, 1/6, and 1/9. For 7 = 1, 1/4, and 
1/9, we considered three different values of Mt. In total, we used 13 different combinations 
for the two colliding protoplanets (see tabled]). For each mass combination, we varied the 
impact angle in the range = 0° — 75° in 15° steps, and the impact velocity in the range 
Vimp = 1-0 — 3.0fesc in 0.2t>esc stcps, whcrc Vesc is the two-body escape velocity defined as 



(8) 
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where Rt and Ri are the radius of the target (larger protoplanet) and the impactor (smaller 
protoplanet), respectively. To precisely determine the transition between merging and hit- 
and-run collisions, we varied fimp with a smaller step size of 0.02 fesc near the transition. In 
total, we performed more than 1000 runs, consisting of 13 (mass combinations) x 6 (angles) 
X 16 (11 runs with 0.2fesc steps and ~ 5 runs with 0.02fesc steps). 

The impact parameters f imp and 6 are defined when the two protoplanets are in contact 
with each other (see Figured]). We assumed that the two protoplanets are mass points, and 
calculated backward the positions of two mass points until their distance was apart at a 
distance of 3 (Ri+Rt). Then, we performed the giant impact simulations over a period of 
10^ sec. 



3. Collision Outcomes 
3.1. Merging and Hit-and-Run Collisions 

The outcomes of the colhsions between the protoplanets are divided into two types : 



merging and hit-and-run collisions. According to the previous studies ( lAsphaud 12010 



Leinhardt et al.ll2010l : lLeinhardt &: Stewartll201ll ). the collision outcomes are subdivided into 



several regimes (e.g., partial accretion, fragmentation and so on). However, we classify col- 
lision outcomes into only two regimes of merging and hit-and-run collisions, because it is 
the most essential to the evolution of protoplanets during the giant impact stage that only 
one big body is left after the giant impact (i.e., merging collision) or two big bodies are left 
(i.e., hit-and-run collision). Additionally, we need to classify collision outcomes as simply as 
possible, in order to incorporate those into A''-body simulation. 

Figure [2] shows snapshots of two typical collisions. Panels (a) to (h) in Figure [2] show the 
time sequence for a relatively low- velocity collision (fimp = l-3fesc) of same-sized protoplanets 
(Mi = Mt = O.IM®) with 6 = 30°. After the first contact, the protoplanets become separated 
(see panel (e)), but remain gravitationally bound. Although some amount of mantle material 
is ejected, almost all parts of the colliding protoplanets finally merge. We refer to this type 
of collision as "a merging collision" . On the other hand, a relatively high- velocity collision 
leads to a completely different result. Panels (i) to (1) show the time sequence of a collision 
with t>imp = l.Sfcsc- The impact angle and protoplanet masses are the same as in panels 
(a) to (h). After the first contact, the protoplanets escape from each other and are no 
longer gravitationally bound. We refer to this type of collision as "a hit-and-run collision" . 
By examining snapshots at t = lO^sec, we could easily classify almost all the simulations 
performed in this study as either merging or hit-and-run collisions. However, in some cases. 



- 8 - 



the protoplanets were separated but still gravitationally bound at t = 10^ sec. Since those 
protoplanets are expected to eventually merge, we classified such a case as a merging collision. 



3.2. Mass of Protoplanet after a Collision 



Here, we consider the mass of the protoplanet after the collision. We briefiy describe the 
method to determine the mass of the gravitationally bound objects from the SPH particle 
data after t = 10^ sec. In the firs t step, we roughly ideri tify clumps of SPH particles using a 



friends-of- friends algorithm (e.g.. iHuchra fc Gellerlll982l ). We then iteratively check whether 



or not any SPH particles not belonging to clumps are gravitationally bound to clumps. 
Finally, we iteratively identify pairs of clumps that are gravitationally bound. Such clumps 
should merge after t = 10^ sec, and we can thus regard them as a single object. We define 
the mass of the largest object as Mi and the second largest object as M2. 

Figure E] shows the mass of the largest object normalized by the total mass, Mi/M^, as 
a function of the impact velocity normalized by the escape velocity, Vimp/'fesc, for collisions 
of same-sized protoplanet s with 6 = 30°. This figure also shows the numerical results of 
Agnor fc Asphaugl ( l2004a( ) for collisions of same-sized protoplanets with O.IM^, which are in 
good agreement with our results for total mass of 0.2Mq). We also performed simulations 
for Mt = O.4M0 and l.OM^. As seen in Figure[31 Mi/M^ does not depend on the total mass, 
when the normahzed impact velocity, fimp/'^^esc, is considered. The insensitivity to the total 
mass holds true to collisions between protoplanets, which has been predicted by lAsphaug 
( 120101 ). If the material properties such as strength dominates over gravity (typically collision 
between smaller b odies with less t han 1 km in radius), the collision outcomes would depend 
on the total mass ( lAsphaugll2010l ). 



3.3. Transition betvi^een Merging and Hit-and-Run Collisions 



Figure [3] indicates that Mi /Mr changes sharply around 1.4fesc- Collisions at impact 
velocities less than this velocity result in almost perfect accretion (i.e., Mi/Mt — 1), and thus 
are classified as merging collisions. For the case of a near head-on collision {6 < 15°), Mi/Mx 
gradually decreases with the increase of fimp- Although such a coUisior i should be clas sified 
as a partial accretion collisi on or a fragmentation collision according to lAsphaugl ( 120101 ) and 
Leinhardt fc StewartI ( 120111 ). we here classify those as a merging collision because only one 
large body remains after the collision. Since the probability of a near head-on collision with 
high velocity is quite low during the giant impact stage, our treatment would not become a 
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serious problem. On the other hand, Figure [3] indicates that the colhsions at impact velocities 
higher than 1.4fesc result in Mi/Mx — 0.5. Some amount of material become stripped from 
the protoplanets, and the protoplanets escape from each other. These collisions are classified 
as hit-and-run collisions. We refer to the impact velocity at the transition between merging 
and hit-and-run collisions as the critical impact velocity, Vcv The normalized critical impact 
velocity, Vcr/vcsc, is not strongly dependent on the total mass, Mt. We obtain fcr/fesc = 
1.39 ± 0.01, 1.37 ± 0.01, and 1.37 ± 0.01 for Mt = 0.2, 0.4 and l.OM®, respectively. 

The critical impact velocity is expected to depend on the impact angle. Figure S] is 
similar to Figure [3l but for the collisions with 6 = 60°. As is the same in the case of ^ = 30°, 
Mi/Mp for 6 = 60° does not depend on the total mass. However, collisions with 6 = 60° 
result in lower fcr/'^esc values than the 9 = 30° case. The calculated values are fcr/^^esc = 
l.lliO.Ol, 1.09±0.01, and 1.09±0.01 for Mt = 0.2, 0.4 and l.OM®, respectively. This result 
implies that collisions at higher impact angles are more like to be hit-and-run collisions. 

The critical impact velocity is also expected to depend on the mass ratio of the proto- 
planets (7). Figure shows the results for 7 = 1/4 (mass ratio of 1:4) and 6 = 30°. It can 
be seen that in this case also, Mi/Mt does not depend on the total mass. This has also 
been verified for 7 = 1/9 (mass ratio of 1:9), although the results are not shown here. As 
shown in Figure [5l Mi/Mt changes sharply around l.Sfesc- The normalized critical impact 
velocities are v^Jvesc = 1.57±0.01, 1.53±0.01, and 1.53±0.01 for Mt = 0.5, 1.0 and 1.5Me, 
respectively, which are larger than for the case of collisions between same-sized protoplanets. 



3.4. Dependence on Particle Number and Initial Internal Energy 

In addition to the impact simulations with 20,000 SPH particles (standard case), we 
performed simulations with 3,000 (low-resolution case), 60,000 and 100,000 particles (high- 
resolution cases) for certain impact parameters in order to check the dependence of con- 
vergence on particle number. Although we found that the critical impact velocity for the 
low-resolution case was slightly different from that for the standard case, the high-resolution 
cases yielded the same results. For example, for 7 = 1, Mt = 0.2M^, and 6 = 30°, we 
obtained fcr/^csc = 1-39 ± 0.01 for both the standard and high-resolution cases. Therefore, 
using 20,000 SPH particles is enough for determining the critical impact velocity. 

In subsection 12. 2[ the initial internal energy of the SPH particles was set to 1.0 x 10^ 
J/kg. To investigate the effect of the initial thermal state of the protoplanets, we prepared 
pre-impact protoplanets with an internal energies of 1.0 x 10^ (cold-state case) and 3.0 x 
10® J/kg (hot-state case). We then performed simulations of collisions between same-sized 
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protoplanets with impact parameters near the transition between merging and hit-and-run 
colhsions. We found that although the results for the very low-impact angle case {6 = 15°) 
showed a very slight dependence on the initial thermal state, for all other cases no such 
dependence was observed. 



4. Merging Criteria for Colliding Protoplanets 



In the previous section, we determined the critical impact velocities for several impact 
parameters. In Figure El we summarize Vcv/vesc for all parameter sets of the mass ratios (7) 
and impact angles (6). For 7 = 1, 1/4, and 1/9, we performed simulations for three sets of 
Mt (see Tabled]). In Figure [6l we plot the average of these three results. 



The critical impact velocities estir nated bvlAgnor fc Asphaud (l2004al ) and lAgnor fc Asphaug 



( l2004bl ) are also plotted in Figure El lAgnor fc Asphaugl ( l2004al ) investigated collisions be- 
tween same-sized protoplanets with masses of O.IM0, and found that Vcr/vcsc = 1-4 — 1.5 
for 9 = 30°, and Vrr/vef^r. = 1-1 ~ 1-2 for 6 = 45° or 60°. Our results are in good agreement 
with those of lAgnor &: Asphaugl (l2004al ). since in our simulations, Vcr/vesc = 1-39 ± 0.01 for 
9 = 30°, 1.19 ± 0.01 for 9 = 45°, and 1.11 ± 0.01 for 9 = 60°. 



Agnor fc Asphaugl ( l2004bl ) investigated collisions be tween different- sized protoplanets 



with mass ratios of 1:2 and 1:10. Based on Figure 17 in lAsphaud ( l2009l ). for a mass ratio 
of 1:2, t;cr/t^csc = 1-4 - 1.5 for 9 = 30°, and 1.1 - 1.2 for 9 = 45° or 60°. In the present 
study, for a mass ratio of 1:2, we obtained fcr/t'csc = 1-43 ± 0.01 for 9 = 30°, 1.19 ± 0.01 
for 9 = 45°, and 1.09 =b . 01 for 9 = 60°, and these values are in good agreement with 
those of lAgnor fc Asphaugl (l2004bl ). Although we did not perform simulations for a mass 
ratio of 1:10, our d ata po ints for 1:9 (7 = 1/9) seem to fall within the range obtained by 
Agnor fc Asphaugl J2004b[ ). 



As shown in Figure El fcr/^^esc increases with decreasing impact angle or mass ratio, 
which means that collisions with low impact angles or low mass ratios tend to be merging 
events. This can be explained in terms of the size of the overlapping volume of the colliding 
protoplanets. Since this is geometrically smaller for higher impact angle, the fraction of ki- 
netic energy converted to thermal energy of protoplanets and kinetic energy of the fragments 
is small, resulting in a hit-and-run collision. In addition, in the case of a small impactor (i.e., 
small 7), most of the volume of the impactor tends to overlap with the target. Therefore, 
the impactor can not easily be ejected, which leads to be a merging collision. 



In the following, we consider a simple physical model in order to express the critical 
impact velocity as a function of the impact angle and mass ratio. For two spheres with 
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radii -Rt and Ri colliding with an impact angle 6 (see Figure [7]), the mass fractions of the 
overlapping volumes for the target and impactor (/3t and f3i, respectively) are geometrically 
given by 

A = ^=(-^j (1-sm^) |---^^_j(l-sm^)j, (9) 



Mll=[ (^)'(l-sin^)2||-i(^)(l-sin^)|, if sin0>l-jg^, 

otherwise, 

(10) 

where and are the masses of the overlapping volumes for the target and impactor, 
respectively. A constant density is assumed for simplicity. 

As illustrated in Figure [TJ we divide the spheres into overlapping parts and non- 
overlapping parts, and consider the momentum exchange between the overlapping parts. 
Using the parameter of the degree of the momentum exchange (a), the post impact veloci- 
ties of the overlapping parts (f "post "^i^post) expressed as 

■^"post = (1 ~ Ol){Vt — ■WcOm) + "^COM) /^^\ 
<;ost = (l-«)(^i-^COM)+^COM, ^ ' 

where fcoM velocity of the center of mass of the overlapping parts, and is written as 

COM j^ov ^ov • ) 

For example, when a = 1, the velocities of the overlapping parts (fj^^post f°post) become 
'^COM owing to the complete momentum exchange. 

The post impact velocities of the entire target and impactor can be derived based on 
the conservation of momentum, and written as 



MiVi^post = (M - Mr)vi + M^-u 



ov „ ,ov 
i,post; 



where ft,post and fj^post are the velocities for the post-impact target and impactor, respec- 
tivery. Combining equations. ( ITT]) - (IT^ gives 



(14) 
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If the relative velocity of post-impact objects is higher than the two-body escape velocity, a 
hit-and-run collision should occur. Therefore, to solve f i,post — ^t,post = "^^esc using = fi — ft, 
the critical impact velocity is obtained as follows: 



^esc 



aPtPi- 



'AMt + AMi. ■ ^^^^ 

The calculated results for a = 0.6 are drawn as gray curves in Figure [H We could roughly 
reproduce the dependence of fcr/^^esc on 9 and 7, but detail features, especially the case for 
the low-impact angle and low-mass ratio, could not be reproduced. For such a collision, the 
role of fragmentation that has not been considered here may become important. It should be 
also noted that a = 0.6 would not be applied to collisions other than giant impacts between 
rocky protoplanets. 

In addition to the above physical model, we tried to mathematically fit the numerical 
data for the critical impact velocity. Of the many possibilities available, we found that the 
following simple formula with five fitting parameters, Ci to C5, was most effective, 

^ = ciFO'^ + c,r + c^e''^ + C4, (16) 

"^csc 

where F = (1 — 7)/(l + 7) = (Mt — Mi)/MT, and Q = \ — sin Q. The fitting parameters are 
ci = 2.43, C2 = —0.0408, C3 = 1.86, C4 = 1.08, and C5 = 5/2. The fitting curves produced 
by equation f|T6|) are shown in Figure [8] as thick curves, and are in excellent agreement with 
numerical results. The value of the fitting parameters derived here is limited to the collision 
between protoplanets. 



5. Merging Probability of Protoplanets 



Almost all previous N-body simulations of terrestrial planet formation during the giant 
impact stage have been performed based on the assumption of perfect accretion. In o r der to 
investigate the statistical properties of fully formed terrestrial planets, iKokubo et al.l (120061 ) 
considered 10 sets of protoplanet initial conditions, and performed 20 runs for each set under 
the assu mption of perfe c t accr etion. Subsequently, to investigate the spin state of the formed 
planets, iKokubo fc Idal (120071 ) performed additional 30 runs each for 7 sets of protoplanet 
initial conditions. Using the formula for the merging criteria (equation [16]) derived in the 
present study, we can now determine whether each of giant impacts was a merging or a 
hit-and-run event. 



Figure [9] shows the normalized impact velocity as a function of impact angle (left panel) 
and mass ratio (right panel) for 635 giant impact events during 50 runs under standard initial 
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conditions (Model 1 in iKokubo &: Idal (120071 )). The symbols denoted by crosses represent 
hit-and-run events, as determined by equation f|T6|l . In fact, 40% of all impact events (256 



out of 63 5) are expected to be hit-a. nd-run collisions. This result is consistent with a previous 
study by lAgnor fc Asphaugi (l2004a( ). who estimated a hit-and-run probability of roughly half. 



6. Summary and Discussion 

During the giant impact stage of terrestrial planet formation in our solar system, a few 
tens of Mars-sized protoplanets collide with each other to form terrestrial planets. Almost 
all previous studies on N-body calculations of the giant impact stage have been based on 
the assumption of perfect accretion. However, recent impact simulations have shown that 
collisions of protoplanets are not always merging events. 

As a first step towards studying the effects of such imperfect accretion on terrestrial 
planet formation, we investigated the merging criteria for a collision of rocky protoplanets. 
Using the SPH method, we performed more than 1000 simulations of giant impacts for various 
parameter sets, such as the mass ratio of colliding protoplanets (7), the total mass of two 
protoplanets (Mt), the impact angle (6*), and the impact velocity (wimp)- We investigated the 
critical impact velocity (fcr) at the transition between merging and hit-and-run collisions. 
We found that the normalized critical impact velocity, fcr/^^esc, depends on 7 and 0, but does 
not depend on Mt. We derived a simple formula for fcr/fcsc as a fun ction of 7 and 6 (see 



equ ation [TBI), and applie d it to the giant impact events considered by iKokubo et al.l (120061 ) 



and lKokubo fc Idal (120071 ). We found that 40% of these events should not be merging events. 



Kokubo fc Gendal (120101 ) was the first to performed N-body simulations of the giant im- 
pact stage taking into account the merging criteria shown in equation (|T6|) . and investigated 
the effects of imperfect accretion on terrestrial planet formation. They found that some 
basic properties such as the final number, mass, orbital elements, and growth timescale of 
planets did not change very much, but the spin angular velocity of the fully formed planets 
was about 30% smaller than that for the perfect accretion model. They also determined 
that 49% of all impact events were hit-and-run collisions, which is also consistent with our 
estimate. 

In this paper, we focused on the merging criteria for protoplanet collisions. In the 
future, we plan to investigate additional collisional phenomena such as mantle stripping and 
ejection of small particles. Using our simulation data for more than 1000 collisional events, 
we can estimate the change in the core-mantle ratio during the giant impact stage. This is 
highly relevant to the formation of Mercury, and the formation probability of such planets 
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with very large cores will be investigated. 

The ejection of small particles during each collision in the giant impact stage may 
influence the orbital evolution of terrestrial planets. The ejected material may have damped 
the eccentricities of the terrestrial planets to their present low values, although it depends 
on the total amount of material ejected during the giant impact stage. 
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Fig. 1. — Geometry of the collision between a larger (target) and smaller (impactor) pro- 
toplanets with radii of -Rt and Ri, respectively. Since a center of mass coordinate system is 
used, the impact velocity, fimp, is given by |fi — ft|, where Vt and Vi are the velocities of the 
larger and smaller protoplanets, respectively. The impact angle is 6. A head-on collision 
corresponds to 6* = 0°. 
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Fig. 2. — Snapshots of two typical giant impacts between equal-mass protoplanets with 
O.IM^. Panels (a) to (h) show the time sequence for a relatively low- velocity collision of 
protoplanets (fimp = l-3fesc) with the impact angle 30°. The colliding protoplanets finally 
merge so that this type of collision is referred to as a merging collision. Panels (i) to (1) 
show the time sequence for a relatively high- velocity collision (fimp = l-5fesc)- The impact 
angle and mass of the protoplanets are the same as in (a) to (h), but the protoplanets do not 
merge. This type of collision is referred to as a hit-and-run collision. The dark gray circle 
at the top-right corner in panels (i) to (1) indicates the size of the initial protoplanets. 
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Fig. 3. — Normalized mass of the largest gravitationally bound object, Mi/M^, as a function 
of normalized impact velocity, fimp/i'csc- Data for 6 = 30° impacts between two equal-mass 
protoplanets with dif ferent total masses are plot ted. Our collision outcomes are very similar 
to those obtained by lAgnor Sz Asphaugl ( l2004al ) who performed simulations for collisions of 
same-sized protoplanets with masses of O.IM^. 
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Fig. 4. — The same as Figure [3] but for 9 = 60°. 
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Fig. 5. — The same as Figure [3] but for collisions of different-sized protoplanets with the 
mass ratio of 1:4 (7 = 1/4). The mass of the second-largest object, M2, is also plotted. 
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Fig. 6. — Critical impact velocities for the cases of the various impact angles and mass 
ratios of protoplanets. Note that (Mt — M-^jM^ is a function of 7, i.e., (1 — 7)/(l + 7). 
Data points represent our numerical results for Q = 15° (filled circles), 30° (filled triangles), 
45° (filled squares), 60° (open circles), a nd 75° (open trian g les). B ars l abeled with AA04a 
and AA 04b are the results obtained by Agnor fc Asphaug ( 2004a ) and Agnor &: Asphau^ 
(l2004bl ). respectively. 
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(a) pre-impact (b) momentum exchange (c) post-impact 



Fig. 7. — Configuration of a target and an impactor for tlie simple physical model, (a) 
Before the impact, the velocities of the target and the impactor are ft and V\, respectively, 
(b) During the impact, the overlapping parts (shaded areas) exchange momentum, and the 
post-impact velocity of the overlapping area is reduced to f °post the target and f °post 
for the impactor. (c) After the impact, velocities of the target and impactor are ft,post and 
^^i,post, respectively. 
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Fig. 8. — Fits to the normalized critical impact velocity data. Data points are our numerical 
results and the same as shown in FigureO Thin and thick curves represent fits using equation 
(ITSll with a = 0.6, and equation ( fT6il . respectively. 
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Fig. 9. — Normalized imp act velocity for 635 giant impacts reported by lKokubo et al.l ( 120061 ) 
and iKokubo Sz Ida! (120071 ) as a function of the impact angle (a) and the mass ratio of proto- 
planets (b). Based on the critical impact velocity (equation [ISj), the giant impact events are 
distinguis hed as hit-and- r un co llisio ns (cross symbo l s) or m erging collisions (circle symbols). 
Although iKokubo et al.l ( 120061 ) and iKokubo &: Idal (120071 ) assumed the perfect accretion of 
protoplanets in their N-body calculations, the present study reveals that 40% of the impact 
events (256 out of 635) are hit-and-run collisions. 
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